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We demonstrate a qubit readout scheme that exploits the Jaynes-Cummings nonlinearity of a 
superconducting cavity coupled to transmon qubits. We find that in the strongly-driven dispersive 
regime of this system, there is the unexpected onset of a high-transmission "bright" state at a critical 
power which depends sensitively on the initial qubit state. A simple and robust measurement 
protocol exploiting this effect achieves a single-shot fidelity of 87% using a conventional sample 
design and experimental setup, and at least 61% fidelity to joint correlations of three qubits. 



Circuit quantum electrodynamics (cQED) is the study 
of the interaction of hght and matter where supercon- 
ducting qubits playing the role of atoms are strongly cou- 
pled to microwave transmission line; resonators acting as 
cavities [1] . This architecture offers advantages over cav- 
ity QED with atomic systems with regard to both cou- 
pling strength and accessibility of strongly driven nonlin- 
ear regimes [2-4]. For example, the high-power response 
of a cavity with a qubit in resonance has recently been 
studied, demonstrating the appearance of additional pho- 
ton number (n) resonances with characteristic ^Jn spac- 
ing, in excellent agreement with theory [4] . In this Letter, 
we investigate the behavior of a strongly driven cQED 
system in the dispersive regime, where four strongly cou- 
pled transmon qubits [5] are far detuned from the cav- 
ity. We find that for increasing drive strength, the cavity 
response initially becomes strongly nonlinear and con- 
tinuously shifts down in frequency. At a critical power, 
it reaches its bare frequency and transmission rapidly 
rises to a "bright state" . Importantly, this critical power 
is strongly dependent on the initial state of the qubit, 
providing for a simple high-fidelity qubit readout mech- 
anism. This scheme is also applicable to the joint mea- 
surement of several qubits simultaneously. 

There are two commonly used cQED readout schemes. 

The first employs the state-dependent dispersive shift of 
the coupled cavity's resonance frequency [Fig. 1 (a)] to in- 
fer the qubit state by measuring transmission [6, 7]. This 
must be done in the low-power, linear response regime of 
the cavity in part because system anharmonicity inhibits 
higher population. The signal-to-noise ratio (SNR) and 
fidelity of this weak measurement is thus limited (to typi- 
cally 40%-60% [7, 8]) by the relatively high noise temper- 
ature {T^^ K, 10 K ~ 20 photons) of conventional com- 
mercial amplifiers and short integration time mandated 
by qubit relaxation (Ti « 1 jis) [9]. A specialized super- 
conducting low-noise amplifier [10-12] could improve the 
SNR, but at a cost of additional complexity. The second 
scheme uses some nonlinearity to project the qubit state 
onto a classically distinguishable system and yield high 
discrimination [13]. The Josephson bifurcation amplifier 
(JBA) is one such example, exploiting cavity bistability 



caused by the Kcrr-DufRng nonlinearity of an additional 
Josephson junction [14-16]. These schemes have, how- 
ever, not yet taken advantage of the nonlinearity inherent 
to the dispersive interaction. The mechanism we report 
here also projects the qubit state onto a classical state 
of the system but does not require additional Josephson 
junctions or any other hardware modifications. It instead 
employs the nonlinearity due to the Jaynes-Cummings 
interaction of the qubit and cavity, in qualitative agree- 
ment with theory [17, 18], effectively using the qubit as 
its own amplifier. 

Our readout scheme exploits the unusual behav- 
ior of the system when driven strongly. For low 
drive strength, the cavity (at /bare = 9.070 GHz) is 
dispersively shifted many linewidths by the strongly 
coupled transmons at frequencies (/i, /2, /s, .A) = 
(6.000, 7.000, 8.000, 12.271) GHz ± 2 MHz (see Supple- 
ment for device details). The cavity also inherits anhar- 
monicity (a, the frequency difference between subsequent 
cavity transitions) from the interaction, which inhibits 
transmission for larger power [Figs, l(a-d)]. For small 
number of excitations n, a is to first order constant [19], 
though higher order terms are often considered as n ap- 
proaches a critical value [20]. This expansion neglects 
the high-power behavior of the system, which can be re- 
covered with a semiclassical model, shown previously for 
the case of a qubit and cavity in resonance [21]. In a 
parallel publication [17], Bishop et al. consider the case 
of a large cavity-qubit detuning scmiclassically and show 
that for increasing n, a decreases monotonically and the 
cavity frequency (/c) shifts toward /bare, in qualitative 
agreement with our observations. The measured cavity 
output power as a function of incident drive power and 
frequency thus smoothly evolves from a linear relation at 
low occupation (n ^ 1), where the shift due to increas- 
ing n is small compared to cavity linewidth, to a non- 
linear relation at intermediate n where a is important 
[Fig. 1(f)]. A dip in transmission likely caused by inter- 
ference between two bistable solutions follows the cavity 
down in frequency with increasing power and disappears 
at /bare- At large drive power, the cavity reaches a near- 
unity transmission "bright state" where /c — > /bare and 




FIG. 1. (Color) Cavity transmission as a function of drive power and frequency, demonstrating tiie "bright" state at high 
incident power, (a) Dispersively shifted cavity response for excited (blue) and ground (red) states of the 8 GHz qubit with 
~ 1 photon mean cavity occupation. We reference this power to dB. For all plots in this figure, the 8 GHz qubit is prepared, 
a measurement tone is pulsed on, and the responding homodyne amplitude is averaged for 400 ns to yield Ah- The mV scale 
used is arbitrary, but consistent to ease comparison. The a;-axis for (a-e) is frequency and covers the same range as for (f-g). 
(b-d) Cavity response for increasing drive power, with data for previous power plotted with dashed lines. Transmission is 
inhibited by the cavity's inherited nonlinearity, limiting dispersive measurement fidelity. Note for example that increasing drive 
10 dB from (b) to (c) increases by only a factor of ~ 2 and complicates the frequency dependence. The emergence of a 
distinct resonance can be seen in (d). (e) High-transmission bright state. At large drive power, the Jaynes-Cummings cavity 
anharmonicity shrinks sufficiently to allow near- unity transmission at /bare. For this power the system only reaches its bright 
state when the qubit is excited due to the asymmetry of the dispersive cavity shift about /bare. This asymmetry is characteristic 
to the transmon qubit [5] but might be possible to simulate for other designs [17]. (f-g) Cavity response (log magnitude) for 
qubit ground and excited states. The cavity continuously evolves from its low power linear behavior through the anharmonic 
bistable region and to the bright state. There are two peaks present in (g) due to qubit relaxation during measurement. The 
symbols {+) denote the optimal power and frequency for qubit readout, where the cavity response for the two qubit states is 
maximally different, (h) Response at /bare (arrows) versus input power, showing a steep jump in transmission corresponding 
to the onset of the bright state at a qubit-state-dependent power. Though transmission state dependence exists elsewhere, the 
behavior here is especially amenable for use as a qubit readout because the difference is large compared to amplifier noise. 



a —7- 0. This power is ^ 50, 000 times larger than what 
induces a one-photon mean population, producing a sig- 
nal large compared to measurement noise. We reference 
the one-photon power to dB. 

This behavior can be used to realize a high-fidelity 
qubit measurement because of a qubit-state-dependent 
bright-state onset power. When driving with sufficiently 
high power at /bare 7 the system will excite to its bright 
state and remain there as long as the drive is applied. 
Crucially, however, the lowest power for which the system 
will excite depends strongly on the qubit state because 
the cavity is dispersively shifted closer to /bare when the 
qubit is excited [Fig. 1(g)]. This smaller detuning reduces 
the Lorentzian filtering effect of the cavity at /bare, low- 
ering the onset power. The cavity response at /bare thus 
undergoes a sharp step for increasing drive at a power 



depending on the qubit state [Fig. 1(h)]. This transition 
does not happen discontinuously; there are powers where 
the system will only excite some of the time. A high- 
discrimination readout is then made by choosing a drive 
power Prncas that maximizes the difference in bright-state 
probability for the qubit ground and excited states, gen- 
erating a large state-dependent transmission difference 
[Fig. 2(a)]. For measurements of the 8 GHz qubit, the 
optimal power is Prncas ~ +47 dB. 

We determine the readout measurement fidelity by 
comparing the ensemble distribution of cavity responses 
having prepared the qubit ground and excited states. 
This is done by integrating ~ 500,000 (in this case) in- 
dividual measurement transients such as those shown in 
Fig. 2(a) for 500 ns for both qubit preparations and plot- 
ting their distribution [Fig. 2(b)]. For the qubit at 8 GHz, 
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FIG. 2. Measurement transients and histograms for single- 
shot measurement of one qubit (at 8 GHz), (a) Three mea- 
surements of single-shot system response Vnit) to driving at 
/bare with a power Pmeas for prepared qubit ground and ex- 
cited states, low-passed to the response time of the cavity 
1/k « 100 ns. They are reproducible and well-distinguished, 
demonstrating that the response is large compared to the 
measurement noise, (b) Histograms and S-curves quanti- 
fying measurement fidelity. An ensemble of single-shot re- 
sponses are integrated for 500 ns and their distribution plot- 
ted. The two histogram peaks (solid lines) are well separated 
with few counts between them. Integrating these yields "S- 
curves" (dashed lines), with their maximal difference indicat- 
ing a single-shot fidelity of 87%. (c) Measurement fidelity 
versus integration time. Integrating Vnit) for 120 ns yields 
80% fidelity, while 240 ns yields 87%. The histogram evolu- 
tion shows that integrating longer further separates the peaks 
but does not increase fidelity. 



the maximal separation between the normalized integral 
of these histograms indicates a measurement fidelity of 
87%. This increases to 91% when pulsing to the second 
transmon excited state just before measurement (data 
not shown; similar to Ref. 16). Fidelity is likely limited 
by qubit relaxation before its state is projected (false 
negatives) and spurious system excitations due to insuf- 
ficiently separated onset powers (false positives). Quan- 
titative predictions of fidelity based on parameters like 
qubit lifetime are unavailable because the fast dynamics 
of the measurement process are not yet well understood. 

Two important considerations for qubit readout are 
measurement time and repetition rate. Single-shot fi- 
delity as a function of integration time is shown in 
Fig. 2(c). After 120 ns the fidelity has risen to 80%, and 
by 240 ns, it reaches its ultimate value of 87%. This time 
is due to the speed with which the system excites to the 
bright state and the SNR. Integrating longer further sep- 
arates the histogram distributions but does not change 
the fidelity because it is limited by failed encoding of the 
qubit state onto the cavity, not noise. The large number 
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FIG. 3. Joint qubit readout. (a) Pulsed cavity re- 
sponse Au when preparing the |111) state. The eight 
(2^) register states of three qubits induce a different 
dispersive cavity shift, each discernible at low power 
due to decay of the |111) state during measurement. 
These frequencies were independently measured to be 

(/|000) , /|100), /|010) , /|001) , ./lllO), /|101), /|011), /|111>) ~ 

(9.145, 9.139, 9.131, 9.112, 9.124, 9.105, 9.097, 9.090) GHz, 
with /bare ~ 9.070 GHz. The most prominent secondary 
cavity position corresponds to the third qubit (at 8 GHz) 
relaxing during measurement, consistent with its Purcell- 
limited [22] Ti being shortest of the three, with T^^ = 1.2 fis, 
yQ2 _ ]^ Q gj^^ rpQs _ Q g rpj^g system excites to its 
bright state at lower power here than seen in Fig. 1 because 
of the smaller initial cavity detuning from /bare. The color 
scale is identical to that in Fig. 1. (b) Histograms for all eight 
basis states, when driving with -1-49 dB at /bare, (c) S-curves 
for joint qubit readout, indicating a minimum of 61% fidelity 
for distinguishing the |000) from any other state. The best 
fidelity is in this case obtained with a ~ 25% probability of 
exciting the system when prepared in |000). 



of measurement photons used might be of concern re- 
garding the creation of quasiparticles in or heating of the 
sample. We have not observed any effect on qubit coher- 
ence when using this scheme, however, suggesting that 
any effects are gone after the 10 fxs repetition time, nor 
does its use affect the base temperature of our ■^He-'*He 
dilution refrigerator (25 mK). An example measurement 
of a Rabi oscillation is shown in Supplementary Fig. SI. 

We have extended this readout scheme to the joint 
readout of three qubits. The physics of multiqubit read- 
out are similar to the solitary case, but generalized for the 
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eight zero-power cavity frequencies possible with three 
qubits. Each basis state induces a different dispersive 
cavity shift from /bare; yielding a hierarchy of powers for 
the onset of the bright state. Readout is done by choos- 
ing a measurement power above the onset power for all 
but the state one is trying to detect. We consider the reg- 
ister of three qubits at frequencies (6, 7, 8) GHz, with the 
fourth qubit unused. For these detunings, Pjoint 

« 49 dB 

is optimal for distinguishing |000) from all other basis 

states (including second excited states). This scheme is 
especially convenient for performing qubit tomography, 
which involves measuring multiqubit correlations [23] and 
is done efficiently with a joint qubit measurement [24- 
26]. Indeed, the creation and detection of tripartite en- 
tanglement was recently demonstrated using this device 
and readout scheme [27]. Cavity transmission for the 
prepared excited states of all three qubits is shown in 
Fig. 3(a). 

The joint fidelity is quantified in a similar manner as 

for one qiibit, but expanded for the eight qubit basis 
states [Fig. 3(b)]. As there is no established definition 
of joint qubit readout fidelity, we propose a conservative 
metric: the maximal separation between the S-curvcs for 
the state one is trying to detect and the least distinguish- 
able state. This standard yields a fidelity of 61% to the 
|000) state for these qubit detunings [Fig. 3(c)], limited 
by the distinguishability of the excited state of the first 
qubit (at 6 GHz). Fidelity could be substantially im- 
proved by cither reducing that qubit detuning (at a cost 
of increased residual qubit-qubit coupling) or by pulsing 
to the second transmon excited state. When an addi- 
tional TT pulse between the first and second excited states 
is included, the joint fidelity is increased to ~ 80% (data 
not shown). This scheme could presumably be scaled to 
a larger number of qubits, subject to concerns of spectral 
crowding and residual qubit coupling. 

This measurement scheme, which we name the Jaynes- 
Cummings readout, provides high-fidelity single-shot 
qubit measurement in cQED while being simple, rela- 
tively insensitive to qubit detuning, and not requiring 
any change in experimental setup or sample design. It 
works by combining the occupation-dependent Jaynes- 
Cummings anharmonicity [17, 18, 28] with the qubit- 
state-dependent dispersive cavity shift to conditionally 
drive the system to a high-transmission bright state, in 
analogy to the fluorescence readout scheme used in ion 
traps [29]. The effect has also delivered similar perfor- 
mance in a Purcell-filtercd [30] low-Q cavity coupled to a 
single transmon, suggesting that it is also robust to cav- 
ity lifetime. It may also be extendable to other cQED 
architectures, not just those using transmons [17]. We 
do not expect it to be quantum nondemolition due to 
dressed dephasing [31] of the qubit from the substantial 
photon population, but this is not necessary for many ap- 
plications. Future work will further optimize the readout 



scheme by changing sample design parameters. A scheme 
may be possible where all qubit states are measured in 
a single shot by sweeping the measurement tone power 
and frequency or by using an intermediate power where 
the transition rate for each basis state is substantially 
different. 
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Device Fabrication and Parameters 

The device is fabricated on a C-plane sapphire (AI2O3) 
substrate. Optical hthography followed by fluorine-based 
reactive ion etching is used to pattern the coplanar 
waveguide (CPW) structure comprising the resonator 
onto sputtered niobium. The resonance frequency of this 
structure, in the absence of the qubit dispersive shift, is 
the so-called "bare" cavity frequency /bare = 9.070 GHz. 
Gaps in the CPW at either end of the chip comprise cou- 
pling capacitors to the environment, giving the cavity a 
Hnewidth n = 2.4 MHz. 

Electron-beam lithography is used to pattern the 
four transmon qubits, which are deposited using a 
double-angle aluminum electron-beam evaporation pro- 
cess. Maximum Josephson energies for qubits Qi to Q4 
are, respectively, EJ^/h = (42,29,47,57) GHz, with 
qubit-cavity coupling .g/27r sa 220 MHz (not including 
the multiplicative level-dependent coupling matrix ele- 
ment; see the Methods section of 1) and charging energy 
Ec/h w 330 MHz for each qubit. 

Details of experimental setup and control are available 
in a companion publication [1]. 

Readout Optimization Procedure 

Setting up the Jaynes-Cummings readout in a cQED 
device is straightforward. The bare cavity frequency can 
be first found by measuring transmission through the de- 
vice as a function of frequency, using a very large drive 
power. It is wise to verify that the measurement ampli- 
fier chain will not be saturated by full transmission of 
this signal ahead of time. Alternatively, if the Tc of the 



cavity material is higher than the Tq of the qubit mate- 
rial, the bare cavity frequency can be measured with low 
power at a temperature between their transition temper- 
atures. One should take care to account for the effects of 
kinetic inductance, however, which can be minimized by 
choosing a temperature as low as possible. 

Once /bare is ascertained, the optimal readout power 
for a given qubit detuning can be found in either of two 
ways. The first and more straightforward way is to pre- 
pare the qubit in either the ground or excited state and 
measure transmitted amplitude as a function of drive 
power. The optimal power can be taken to be that which 
max;imizes the difference between the transmission for the 
ground or excited qubit states. In the case of multi-qubit 
readout, each basis state (or just the ones likely to limit 
fidelity) should be prepared and measured. This method 
has the advantage of being fast and easily automated, 
though it does not necessarily optimize measurement fi- 
delity per se, but rather measurement contrast. 

The second way to find the optimal readout power is to 
directly measure readout fidelity as a function of power. 
At each power, an ensemble of single-shot measurement 
transients for both prepared states should be integrated, 
and their fidelities computed as described in the main 
text (see Fig. 2). This method is unlikely to give a 
substantially different optimal drive power than the first 
method and is substantially more involved, but is tech- 
nically the correct way of optimizing the conventional 
definition of measurement fidelity. 



[1] L. DiCarlo et al, Nature 467, 574 (2010). 
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FIG. SI. Quantum measurement following a Rabi oscillation to demonstrate the Jaynes-Cummings readout scheme. The 
qubit is prepared in the superposition state = cos(&/2)|0) + sin(6'/2)|l) and repeatedly measured as a function of 6. In 
this case, 6 is chosen by scaling the amplitude of a gaussian-shaped microwave pulse at the transition frequency of the qubit. 
(a) The distribution of integrated values are shown for 6 = Q, 9 = 7r/2, and 6 = -k. The distribution is always bimodal with 
very few counts between the two peaks for all 6, corresponding to the system either staying in its dark state or becoming 
bright almost immediately, demonstrating the single-shot nature of the readout, (b) Measurement response distribution for 
preparation angles between — —2tt and 6 = 27r. The y-axis is centered on and scaled to the standard deviation (a) of the 
ground state distribution, showing a separation between peaks of about 23a for this integration time of 500 ns. The average 
value is plotted (white dots), showing the canonical sinusoidal Rabi oscillation. 



